A scheme for direct observation of entanglement for Gaussian continuous variables 
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We suggest an experimentally realizable scheme to test entanglement of a mixed Gaussian contin- 
uous variable state. We find that the entanglement condition is simplified for the family of Gaussian 
states which are relevant to experimental realization. The entanglement condition is then shown 
to be directly related to joint homodyne measurements. We show how robust the proposed test of 
entanglement is against imperfect detection efficiency 
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In the current development of quantum information 
processing, entanglement, which is quantum correlation 
between particles, plays a critical role. Recently, the 
generation, manipulation and application of continuous- 
variable entangled states have been extensively studied 
10, H, H Q mainly for Gaussian states. It is because the 
Gaussian continuous variables are extremely similar to 
qubit systems in theoretical treatment and they are the 
continuous-variable states that can be realized in a lab- 
oratory. 

A state is said to be entangled when it is not separable. 
A Gaussian continuous-variable state is separable if and 
only if the partial transpose of its density matrix is non- 
negative Q . This condition is equivalent to the possibil- 
ity of imposing a positive well-defined P function to the 
state after local unitary operations 0]. In this paper, 
we show that the entanglement condition is simplified for 
a family of Gaussian states which arc more relevant to 
experimental realization. We then suggest a scheme to 
test entanglement using joint homodyne measurements. 

Einstein, Podolsky and Rosen (EPR) paradox is closely 
related to entanglement. There have been discussions 
on the EPR argument for continuous variables based on 
Bell's inequality ||. On the other hand, Reid and 
Drummond | [lO| were concerned with the demonstration 
of the EPR paradox itself as distinct from Bell's inequal- 
ities. They derived a noise level which is violated by the 
EPR paradox. Silberhorn et al. proved the gener- 
ation of an EPR entangled state by the indirect mea- 
sure of the sub- vacuum noise in quadrature correlations. 
Even though the relation of EPR paradox to entangle- 
ment is still to be ravelled, it has been generally accepted 
that there is a subtle difference. We show that Reid and 
Drummond's EPR criterion is only a sufficient condition 
for entanglement. 

A two-mode Gaussian state is represented by a Gaus- 
sian Weyl characteristic function 



C(x) = exp(-ixVx'^) 



(1) 



are bosonic operators. The quadrature matrix can in 
fact be written using 2 block matrices Li and L2 for 
local quadrature variables and C and its transpose C"^ 
representing inter-mode correlation. 



V = 



(2) 



Lemma 1. — If the block matrices Li = L2 and C 
are diagonal, i.e., the quadrature matrix of a Gaussian 
continuous- variable state has the following form: 



fni ci 0' 

712 C2 

ci ni 

C2 n2j 



(3) 



where ni or 712 may be smaller than the vacuum limit 1. 
The state is separable if and only if 



5i52 > 1 



(4) 



where 5i = Ui — \ci\ for i — 1,2. 

Proof — By local unitary squeezing operations, the 
matrix (j^) is transformed into 



Vi = 



fn c 0\ 

n c' 

c n 

Vo c' nj 



(5) 



where n — ^Jri\ri2, c — cx^Jn^Jni and c' — ci\/n\ln2- 
The factor n is directly related to the uncertainty prin- 
ciple to satisfy 71 > 1. For the state with the quadrature 
matrix (H), Simon's separability criterion reads 



^)(n^-c'2)>2n2-K2|cc'|-l- 



(6) 



Define the average and the difference of the correlation 
factors, |c| and |c'|: Ca = (|c| -I- |c'|)/2 and Cd — (|c| — 
|c'|)/2. Using the new parameters and c^, Simon's 
criterion can be written as 

[(7. - caf - (1 + c^Mji + c^f - (1 + cl)\ > (7) 



where x = 52,^2), is the transpose matrix of ^here [{n+Ca)'^-{l+cl)] is positive unless n = 1 and c = 



X and V is the quadrature matrix defined as 2{{xi, Xj}) 
The quadrature variables qi and pi for mode i{= 1,2) 
correspond respectively to quadrature operators qi = 
{hi + a[)/\/2 and pi = i{ai — a,|)/\/2 where and aj 



c' = when it becomes zero. The separability condition 
is satisfied if and only if 

(n - Caf - (l-f C^) > ^ (77 - |c|)(77 - |c'|) > 1. (8) 
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With use of the definitions of n, c and c' for the inequahty 
in the right-hand side of the arrow we obtain the separa- 
biHty condition in Eq. (|^). 

One can also use Duan et a/.'s separability criterion |^ 
even though one has to be careful because all the diagonal 
elements in their standard form II have to be larger than 
1 for their test of separability. 

Now, recalling the definition of the quadrature matrix, 
the separability criterion can be written as 

{{q!) + {ql) - 2|(9i<Z2)|)((p?) + {pD 2|(pip2)|) 
> {{ql) + {ql) - 2|(9i<Z2)|)o((p?) + {pD 2|(piP2)|)o 

(9) 

where {X)q denotes mean value of X for the vacuum. 
The right-hand side of the inequality is 1, which can be 
easily seen by substituting n = and s = into (|l^) . We 
have found that the entanglement of a Gaussian field in 
the form (H) can be tested by comparing the quadrature 
correlation of the field with that of the vacuum as shown 
in Eq. (|). 

What is the relevance of the quadrature matrix given 
in (H) to our study of entanglement for continuous vari- 
ables? Among many possible ways to produce Gaussian 
entangled states, two methods are experimentally more 
relevant: One is to use a non-degenerated parametric am- 
plifier (NOPA) to produce a two-mode squeezed state 
and the other is to use a beam splitter as an entangler 
§0. 

Let us consider the entanglement of the output field 
from a beam splitter when two independent Gaussian 
fields are incident on it. The logical definition of the 
density operator ps for a single-mode continuous- variable 
Gaussian state is u% 



pg = J\f exp{-ja'^a la^^ /*a^) 



(10) 



where A/" is a normalization factor (throughout the pa- 
per the normalization factor is generally denoted by M 
even though the detailed forms may differ.). When I = 
and j — huj/kT, where k is the Boltzmann constant, the 
density operator ( |lo|) represents a thermal state of tem- 
perature T. It is straightforward to show that the general 
Gaussian state ( |lo| ) may be transformed into a thermal 
state by the unitary single-mode squeezing operator [ p^ , 
S{Q = exp(|C*a^ — ^C^^^) with the complex squeezing 
parameter ( = se^'^: 



S\C)PsS{C) ^J\fexp{~2^yn^ - \m\^aa^ 



(11) 



for s = {\l\ - j)/4(j + l^l) and ^ = Arg{l). 

Consider that the input field described by the opera- 
tor Si is superposed on the other input field with op- 
erator 02 by a lossless symmetric beam splitter, with 
amplitude reflection and transmission coefficients r and 
t. The output-field annihilation operators are given by 
ci = BaiB'' and £2 = Ba2B'^ where the beam splitter 



operator is 



B 



exp 



-(a|a2e 
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aia^e"*'^) 



(12) 



with the amplitude reflection and transmission coeffi- 
cients t = cos(0/2) and r = sin(6'/2). The phase differ- 
ence 4> between the reflected and transmitted fields can 
be adjusted by putting a phase shifter. 

When the two input fields are Gaussian states, the 
output state from a beam splitter is 



BS, (Ci )52 (C2 )/5f 5] (C2 ) SI (Ci )B 



(13) 



where p*^^ is the density operator for a thermal state and 
the relation (|ri|) has been used. Without losing gener- 
ality, we take the input squeezing parameter to be real 
while keeping </> of the beam splitter variable. 

It has been found by us ||ll| that two squeezed states 
may be maximally entangled when the beam splitter has 
9 = Ti/A and (j) — Q in which case 

BSi{s^)S2{~S2)B^ 

where S'i2(C) = exp(— Caia2 -t- C,*a\a\) is the two-mode 
squeezing operator [ p^ . Throughout the paper, si and 
S2 are assumed to be positive. The two-mode output field 
is represented by the quadrature matrix in the form of 
with its elements 

m = |e-- + |e-= ; .2 = fe^^ + |e-- 

ci = |e---|e-= ; C2 = |e-^ - fe^- 

(15) 

where hi = 2hi -\- 1 with the mean photon number hi of 
the thermal state. Substituting these into the separabil- 
ity criterion (^) we find that the output field is separable 
if and only if 



{hie-"''){h2e-''^) > 1 



(16) 



n^e"'^^* is the quadrature variance of the input field, 
which represents the sub-vacuum noise level when it is 
smaller than 1. 

For the generation of a squeezed state using a NOPA, 
even though it is desirable to produce a two-mode 
squeezed vacuum, it may well be the case that a two- 
mode squeezed thermal field is produced or a two-mode 
squeezed vacuum is produced but decohered in a ther- 
mal bath. Let us thus consider an two-mode squeezed 
thermal field represented by its density operator: 



/5?2(C)=^12(C)/5fpf^l2(C) 



(17) 



By the local rotation operator Ri{(p) = exp(i(/3a|ai/2), 
the two-mode squeezed thermal field can be transformed 
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into pl2{^)- ^ local unitary operation does not change 
the nature of entanglement and the local unitary rota- 
tion is easily realized by a phase shifter. Thus we do 
not lose generality by studying how to test the entan- 
glement of Pi2{s)- We assume that the temperatures of 
thermal fields represented by and /5|'* are the same: 
ni = n2 = n. In this case, the quadrature matrix of 
the squeezed thermal field Pi2{s) has the form (^) with 
matrix elements: 



Entangler 




50:50 Beam splitter 



ni = n2=n,cosh2s 
Cl — — C2 — — nsinh2s. 



(18) 



On the other hand, if a squeezed vacuum is produced 
by a NOPA but decohered in the thermal environment 
with the mean number of thermal photons n, the deco- 
hered state is still represented by the form (||) but with 
its elements |^ 



ni = 712 

Cl -C2 



= e cosh2s -I- n(l — e~ 



-It 



) 



: e sinh2s. 



(19) 



where 7 is the coupling of the field with the environment. 
The influence from the environment to each mode has 
been assumed same. 

We have seen that most of the two-mode Gaussian 
states represented by the quadrature matrices in the form 
(^ are relevant to current experimental techniques. The 
quadrature variables are measured by setting a balanced 
homodyne detector |13, 1^, which is a well-known device 
to detect phase-dependent properties of an optical field, 
at each mode of the two-mode field. The operational 
representation of the balanced homodyne detector is 



Ohd = qcosx- ipsinx 



(20) 



where x depends on the local-oscillator phase. Because 
the phase of the local oscillator is not absolute, we must 
find the local-oscillator phase which gives the off-diagonal 
terms of Li, L2 and C vanish. It is straightforward to 
measure the off-diagonal terms of C by joint homodyne 
measurement as V14 = 2{qip2) and V23 = 2{q2Pi). How- 
ever, measuring the off-diagonal terms of local matrices 
are troublesome as it involves the joint measurement of 
two quadrature variables for a single mode. 

To measure the off-diagonal elements of Li ||l^ , we put 
a 50:50 beam splitter which splits the field in mode 1 as 
schematically shown in Fig. |l|. Using the beam splitter 
operator ( p^ ) for the 50:50 beam splitter, we find that 
the field for three modes 1', 2 and 3 is still Gaussian and 
its quadrature matrix is written as 



, Li + 1 V2C -Li 

-Li + 1 -V2C Li + 1 




(21) 



where the unit matrix 1 is due to the vacuum injected 
into the unused port of the beam splitter. Now the off- 
diagonal elements of Li can be measured by inter-mode 



FIG. 1: Configuration of the setup to test entanglement. The 
dotted devices measure off-diagonal terms of the local quadra- 
ture matrices Li and L2. The boxed device denoted by HD 
are homodyne detectors. The numbers refer to the modes. 



correlation between modes 1' and 3: The mean value of 
the joint measurement (qvPs) — —¥12/^- Similarly other 
off-diagonal terms of the local quadrature matrices can 
be obtained. It is true that if the entangled field is one of 
the fields described in this paper, choosing the phase of 
the local oscillator to vanish the off-diagonal terms V14 
and V23 , all the local off-diagonal terms should vanish but 
to make it sure the above supplementary measurements 
can be performed. 

In fact, as we have shown, we know how to find all the 
matrix elements of the quadrature matrix for a Gaussian 
field so that it is possible to test entanglement not only 
for the fields in the form (^) but also for any Gaussian 
field if the detection efficiency is unity. In this case we 
have to use Simon's separability criterion Q. 

A homodyne detector is composed of two photodetec- 
tors. Inefficient photodetectors introduce noise to each 
mode and reduces the quantum correlation between two 
modes. The detection efficiency may thus determine the 
feasibility of the proposed scheme. If the efficiencies of 
the photodetectors are same, homodyne measurement 
by imperfect detectors is equivalent to homodyne mea- 
surement by perfect detectors following a beam splitter, 
one input port of which is fed by the field to be mea- 
sured and the other by the vacuum |l^. The efficiency 
77 of the homodyne measurement determines the trans- 
mission coefficient ^Jrj of the beam splitter. In fact the 
fictitious beam splitter affects the testing field as though 
it is decohered in the vacuum reservoir. The detection 
efficiency, assumed the same for the both homodyne de- 
tectors, effectively changes the quadrature matrix from 
V to V' = 77V -|- (1 — 77)1. This is what is measured by 
imperfect homodyne detectors. If V is in the form Vq, 
the variance matrix Vq takes the same form as Vq but 
with modified matrix elements n[ = rjUi -I- (1 — 77) and 
|c-| = 77|ci| for each i. 

Consider the effect of the detection efficiency on the 
inseparability of the testing fields. Substituting 77'^ and 
c[ into Eq. (Q) for inefficient detection, we find a state to 
be entangled when 

K - Ic'il)(n2 - 141) - V\h52 - <5i - <52 + 1) 

-1-77(^1 +<52- 2) + 1< 1. (22) 
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FIG. 2: Gaussian states with the symmetric quadrature ma- 
trix Vo in Eq. on the space oi Si = ni — |ci| and 
$2 = n2 — \c2\. Separable states are denoted by S and en- 
tangled states by E and E' with the boundary oi 5iS2 = 1. 
Further the two regions E and E' are separated by the line 
5i+52 = 2. The entanglement imposed on a state of the region 
E' may fail the entanglement test due to inefficient detection. 
However, the entangled state of the region E is always found 
to be entangled regardless of the detection efficiency. 



Rearranging this equation, we can easily find that when 
the original testing field is characterized by (5i + ^2 < 
2, it is always found to be entangled regardless of the 
detection efficiency unless the efficiency is zero. With 
use of (^8|) we find that a two-mode squeezed thermal 
field is extremely robust to the detection efficiency as an 
entangled two-mode squeezed state always passes the test 
of entanglement regardless of detection efficiency. Fig. ^ 
presents the sets of Gaussian states on the space of 6i 
and S2 where separable states with 61S2 > 1 are denoted 
by S and entangled states with 6162 < 1 by E' -I- E'. All 
entangled states in the region E with the condition, 61 -\- 
62 < 2, will violate the inequality (Qj of the quadrature 
correlation unless the detection efficiency is zero while 
some entangled states in E' fail the test of entanglement. 

It is also important to be able to tell a given state is 
pure. A degree of purity can be defined as P = Trp^. 
When P — 1 the state is pure. The purity of a two-mode 
Gaussian state can be performed using ideal homodyne 



detectors. For a Gaussian state of the qu adrature matrix 
the degree of purity P = l/Vdet Vq. Noting that a 
vacuum is pure, the inequality of P < 1 can be written 
as 



2\2 



> 



){{piY - {P1P2Y) 

- {qiq2?)J{p\? ^ {pW2?),.{2i) 



The equality implies that the testing fields are in the pure 
state (with the minimum uncertainty). For inefficient de- 
tectors, once their efficiencies are known as 77, the purity 
may be deduced from the experimental data. 

Reid and Drummond derived the inequality for the 
quantum correlation between two mode fields along the 
line with the Einstein, Podolsky and Rosen argument 
1^, |T^. They introduced the uncertainty Vi (V2) between 
the observable qi {jpi) in one mode and q'l {p'l) inferred 
from the observation of the other mode. Quantum cor- 
relation may lead the product of the uncertainties to be 
less than the vacuum limit, resulting in the inequality of 
V1V2 < 1. In our notation this inequality can be written 



as 



[ni - |ci|)(n2 - IC2I) < 



nin2 



(ni + |ci|)(n2 + \C2\) 



(24) 



Note that the right hand side of the inequality is always 
less than unity. Thus, the inequality (|^) is sufficient to 
satisfy our inseparable condition (see (W). However, the 
converse statement does not hold in general. 

We have proposed an experimental scheme to test the 
entanglement of a Gaussian field for the first time to 
the best of our knowledge. Our scheme is based on the 
inseparability criterion. The present scheme consists of 
balanced homodyne detectors which are well-established 
experimental tools to study quantum optics. We show 
that the entanglement imposed in a Gaussian state is 
measurable persistently against the detection efficiency 
if each mode has the roughly symmetric uncertainty over 
the phase space. 
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